

International Journal of Trend in Scientific 
Research and Development (IJTSRD) 



International Open Access Journal 
ISSN No: 2456 - 6470 | www.ijtsrd.com | Volume - 2 | Issue - 2 


♦ 

♦ 


Thermal Stress Analysis of Homogeneous Orthotropic Plate with 
New Thermal Profiles Using First Order Shear Deformation Theory 


Mr. A. K. Kanase 

Associate Professor, Department of Civil 
Engineering, BSCOER, Narhe, Pune, India 


ABSTRACT 

Laminated composite plates are extensively used in 
Civil, Mechanical and Aerospace industries. A 
composite material is formed by combining two or 
more materials, so that the properties of the 
composites are different from those of the individual 
constituent’s materials. Laminated composites formed 
by fibres are bound together using binding material 
called ‘matrix’. The reinforcing fibers used in fibrous 
composites are having high modulus and ultimate 
strength. The matrix which is used to bind the fibers 
together should be able to transfer the loads to the 
fibers, as well as be able to protect the fibers from 
damage. Analytical models in single layered 
laminated orthotropic plate gives shear deformation 
and transverse normal thermal strains is validated for 
the thermal stress analysis subjected to gradient 
thermal profile across the thickness of laminate. First 
Order Shear Deformation Theory (FOST) is used for 
calculating various quantities and give good 
agreement with gradient thermal profile for thin 
plates. Based on these new results are prepared for 
additional three thermal profiles contains linear as 
well as parabolic path across thickness of plate with 
different boundary conditions. 
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I. INTRODUCTION 

Advancement of the technology of laminated 
materials, it is possible to use these materials in high 
temperature situations now days. Composites have no 
yield-limit, unlike metals and have a variety of failure 
modes, such as fiber failure, matrix cracking, 
interfiber failure and delamination, which give rise to 
a damage growing in service. Composite plates are 
subjected to significant temperature stresses due to 
different thermal properties of the adjacent laminas 
and therefore accurate predictions of thermally 
induced deformations and stresses represent a major 
concern in design of conventional structures. Nature 
of composite plates can be characterized by a complex 
3D state of stress. Various theories and models are 
reported for the thermal stress analyses of the 
laminates. In this article First Order Shear 
Deformation Theory is used for calculating the 
results. Up to todays research scenario most of the 
researchers assumed linear gradient or constant 
thermal profiles along the thickness of plates and 
results for various quantities are validated for these 
thermal profiles. Single layered orthotropic square 
plate is assumed in this paper for preparation of new 
results. Material properties are shown in Table 1. The 
main objective of this paper is to suggest the linear 
and parabolic thermal profiles across thickness of the 
plate with different boundary conditions which is 
shown in Table 2. Thermal profiles variation through 
thickness is shown in Fig.l. Analytical solutions 
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along with additional new results in normalization Table 1 Material Properties for square (a=b) single 
form are prepared. layer orthotropic plate 


First Order Shear Deformation Theory (FOST) 
considered in Reissner, (1945) and Mindlin (1951) for 
solutions to include the temperature effects on 
laminates. Rower et al.,(2001) removed deficiences in 
FOST incorporating third and fifth order displacement 
approximations through the plate thickness. 3D 
elasticity solution can estimate the correct results of 
the thermally induced quantities like displacements 
and stresses. Kant and Swaminathan (2002) 
incorporated formulations through there article 
Analytical solutions for static analysis of laminated 
composite and sandwich plates based on a higher 
order refined theory and suggested First Order Shear 
Deformation Theory as a special case with its 
importance. Kant et al. (2008) explored semi- 
analytical solution for constant and linear temperature 
variation through the thickness of a laminate for 
composites and sandwiches. Kant and Shiyekar 
(2013) developed higher order theory for composite 
laminates subjected to thermal gradient by using 
(HOSNT12) model. 



Material 

Properties 

Ei =150 GPa, Ei = lOGPa, V12 = 0.3, 

1/21 = 0.02 

? 

Gn = G 13 = 5 GPa G23 = 3.378 GPa 
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ax = 0.139E'-6/k a y = 9E-6/k 
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Reference 

T.Kant et al. [5] 

Maximum normalized stress coordinates; a x and a y 
(a/2, b/2,h/2); t (0, 0,h/2). For analysis positive 

coordinates considered. 

Maximum displacement coordinates, u : (0, b/2,h/2); 
v : (a/2, 0, h/2); w: (a/2, b/2,0) 

E\ and Ei : Modulus of Elasticity along laminate 
direction (x-axis) and transverse direction (y-axis). 

V12 andv 2 i : poisons ratio. G 12 , G 13 , G 23 : Modulus of 
rigidity. 

ax,a y : Thermal coefficients along x and y axes. S= 
a/h: aspect ratio. 

a, b, h: Dimensions of plate along x, y and z 
directions. 


Fig. 1 Thermal loading variation of thermal profiles 1 
to 4 across the thickness of the laminate. 
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Table 2 Proposed thermal load profiles along 
thickness of plates 



2 Formulations 
A. Displacement model 

A simply supported single layer orthotropic laminated 
plate is presented along with analytical solution. The 
geometry of the laminate is such that the side ‘a’ is 
along V axis and side 7?’ is on ‘y’ axis. The thickness 
of the laminate is denoted by ‘/P and is coinciding 
with ‘z’ axis. The reference mid-plane of the laminate 
is at h/2 from top or bottom surface of the laminate as 
shown in the Fig. 2. Lamina axes and reference axes 
are coincides with each other. 

The fiber direction of the single layered lamina is 
coinciding with A’ axis of the laminate. Fig.2 also 
illustrates the mid-plane positive set of displacements 
along x-y-z axes. 


B. Constitutive Relationship- 

Stress-strain relationship for orthotropic laminate 
under linear thermal loading can be written as - 
Stress-strain relation for 0° layer . 

Lv>L,, L,) are the stresses and 
{s x ,£y,Y xy , Y yz , y X2 ) are the strains with respect to 

laminate coordinate system (x-y-z). 

dv 

& 

du dw dv dw 

Yxz = — +—,Yyz = 7T + A _ - 
dz dx dz dy 


Where = 


du 

dx 


£ y=- 


_ dw _du dv 
Sz ~~dz’ /xy ~~dy + ~dx 
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III Discussion 

First Order Shear Deformation Theory (FOST) with 
Semi Analytical Model [SAM] presented by Kant et 
al. (2008). Gradient thermal profile 1 is compared 
with results of Kant et al. (2008). From the Table 3 it 
is observed % error for transerse displacement w is 
21.49 for aspect ratio 4 and which decreases to 0.2916 
for aspect ratio 20. For inplane stress ( a x , a ) % error 

decreases from (-3.15,51.31 ) for aspect ratio 4 and 
which decreases to (-2.16,0.12) for aspect ratio 20. 49 
for aspect ratio 4 and which decreases to 0.2916 for 
aspect ratio 20. For shear stresses ( t , r xz , r yz ) % 

^irfi^de&gases from (62.35,37.66,47.61) for aspect 
ratio 4 and which decreases to (3.94,-0.16,0.29) for 
aspect ratio 20. So in all the quantities % error 
decreses as aspect ratio increses from 4,10 and 20. It 
means FOST gives better validation if plate converts 
in to thiner plate. 




Table 3 Maximum normalized displacements () and Stresses (,) for single layered laminated square plate 
simply supported on all edges subjected to gradient thermal loading (TP1) (Aspect ratios 4,10,20) 


s= 

a/h 

SO 

u 

V 

w 



Gy 

Gz 

G z 

4 

TPi 

-0.0811 

-0.2865 

10.9618 

[21.49] 

1.9528 

[-3.15] 

-3.1077 

[51.31] 

-0.5774 

[62.35] 

1.2719 

[37.66] 

-1.2718 

[47.61] 

Kant. [5] 



9.0226 

2.0164 

-2.0538 

-0.3556 

0.9239 

-0.8616 

10 

TPi 

-0.0162 

-0.0228 

1.2226 

[12.93] 

0.4444 

[-7.69] 

-0.5670 

[0.56] 

-0.0613 

[-1.56] 

0.2542 

[-0.89] 

-0.2542 

[0.99] 

Kant.[5]- 



1.4042 

0.4845 

-0.5638 

-0.0638 

0.2565 

-0.2517 

20 

TPi 

-0.0042 

-0.0046 

0.28 

[4.32] 

0.1172 

[-2.16] 

-0.1450 

[0.12] 

-0.0139 

[-3.94] 

0.0659 

[-0.16] 

-0.0659 

[0.29] 

Kant. [5]- 



0.2916 

0.1198 

-0.1448 

-0.0141 

0.066 

-0.0657 


% error= [Kant.(2008)-Present FOST]xl00/Kant.(2008) 
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[ Qa ] are transformed elastic constants or stiffness 


matrix and defined as per the following. 


Q =_£ * 
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a x , a v = Coefficient of thermal expansion in x, y 
direction. (3) 

Compliance matrix involves engineering properties 
namely two extensional moduli (£ ( , E 2 ), two 

Poisson’s ratios ( v n , v 21 ) and three shear moduli 
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Table 4 Maximum normalized displacements ( u , v, 1 OOw) and Stresses (cr x ,<j y , z xy , z xz , z yz ) for single layered 

laminated square plate simply supported on all edges subjected to thermal loading TP 2 , TP 3 and TP 4 . 

(Aspect ratios 4,10,20) 


A/H 

SO 

u 

V 

W 

<7* 



T xz 


4 

tp 2 

-0.040 

-0.575 

5.217 

0.298 

-2.233 

-0.967 

0.635 

-0.635 

tp 3 

-0.036 

-0.345 

4.696 

-0.11 

-2.964 

-0.599 

0.398 

-0.901 

tp 4 

-0.040 

-0.431 

5.218 

0.163 

-2.688 

-0.741 

0.527 

-0.842 

10 

tp 2 

-0.008 

-0.080 

0.601 

0.114 

-0.392 

-0.139 

0.127 

-0.127 

tp 3 

-0.007 

-0.044 

0.541 

0.042 

-0.505 

-0.082 

0.087 

-0.167 

tp 4 

-0.008 

-0.057 

0.601 

0.092 

-0.464 

-0.103 

0.110 

-0.160 

20 

tp 2 

-0.002 

-0.019 

0.139 

0.031 

-0.099 

-0.034 

0.0329 

-0.0329 

tp 3 

-0.002 

-0.011 

0.125 

0.013 

-0.128 

-0.020 

0.023 

-0.0428 

tp 4 

-0.002 

-0.013 

0.139 

0.026 

-0.118 

-0.025 

0.0286 

-0.041 


IV Conclusions 

First Order Shear Deformation Theory (FOST) with 
five degrees of freedom in primary displacement field 
is presented for thermal stress analysis of single layer 
homogeneous orthotropic laminated square plate. 
Thermal loading with doubly sinusoidal variation 
along x and y directions and linear or gradient thermal 
field (Thermal Profill) along thickness direction is 
considered over the cross-ply laminates for validation 
purpose. The present FOST is assessed by comparison 
with Semi Analytical Model [SAM], FOST provide 
good results in case of very thin laminates with aspect 
ratio (S) increases from 4 to 20. Thermal profde 2 
which is gradient thermal profile for same plate only 
boundary conditions of thermal load amplitudes are 
changed as shown in Table 2. So assumption is 
thermal profde 2 having base results for comparison 
of thermal profde 3 and 4. Results canot be compaired 
with each other because all thermal profdes 2,3 and 4 
are totally new therefore new results are presented in 
Table 4 for these additional thermal profdes. Hence 
finaly we predict these additinal thermal loadind path 
along thickness in single layer orthotropic plate 
design. 
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